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Control  of  Nonlinear  Distributed  Paxcimeter  Systems 
with  Application  to  Flow  Control 


The  goal  of  this  project  was  to  put  the  intuitive  idea  of  gain-scheduling  on  a 
rigorous  foundation  for  a  cla^s  of  nonlinear,  distributed-parameter  systems.  This 
involved  a  study  of  the  existence  and  characterization  of  the  ideal,  infinite-dimensional, 
feedback  control.  Since  in  most  applications  the  feedback  function  cannot  be  computed 
in  closed  form  it  was  necessary  to  study  the  convergence  of  approximate  feedback 
functions,  based  on  increasingly  higher  order  finite-dimension«il  approximations  of  the 
system,  to  the  ideal  function.  Finally,  our  results  were  applied  to  Burgers’  Equation, 
which  can  be  viewed  as  a  low-order  approximation  to  a  wide  variety  of  physical 
phenomena,  including  viscous  compressible  flow. 

The  systems  consider  jd  in  this  study  were  modeled  by  an  abstract  differential 
equation  of  the  form 

X  =  F{x,  u) 

where  x  €  9G,  u  €  'll,  the  state  space  9G  and  the  input  space  TL  are  Banach  spaces  with 
'll  finite  dimensional.  In  many  cases  of  interest,  F  will  be  an  unboimded,  nonlinear, 
differential  operator.  For  the  purpose  of  computing  equilibrium  points  and  linearizing 
this  operator  we  are  led  to  restrict  F  to  a  dense  subspace  of  96  with  a  norm  that  makes 
F  differentiable.  By  considering  the  set  of  equilibrium  points 

{(x»,uV))l«e.x} 

we  can  compute  the  linearizations  about  these  points  to  be 

X  =  DjF(x®(a),u®(a))x -I- 

For  each  fixed  a,  we  can  design  a  linecir  feedback  of  the  form 

u  —  G{a)x  -I-  H{a)w 

which  yields  the  desired  closed  loop  linear  system 

X  =  [£>iF(xV), «'’{«))  +  D.,F(x‘'(a),u‘'(a))G(a)]x  -h  D2F{x^{a%u°{a))H{a)w. 

The  objective  of  geiin  scheduling  is  to  derive  a  nonlinear  feedback 

u  =  S{Xyw) 

such  that  the  linearizations  of  the  closed-loop  nonlinear  system  are  equal  to  the  desired 
closed-loop  linear  systems.  This  should  produce  desirable  performance  in  a 
neighborhood  of  the  equilibrium  set.  For  scheduling  to  be  possible  when  the  state 
feedback  G  is  arbitrary,  the  control  transformation  H  must  belong  to  a  family  of 
operators  that  depend  on  G.  We  have  shown  that  a  feedback  of  the  form 


5(a:,  iw)  =  u°(P2{x,  ti;))  +  K{P2{x,  u?))  •  ((x,  w)  -  P(x,  w% 
where  P  =  [Pj  P2]:  95  x ‘U-+9G  x 'll  is  a  projection  from  the  state-control  space  to  the 
equilibrium  set  and  K  ~[G  H],  will  produce  the  desired  result. 

Since  we  cannot  generally  compute  5(x,to),  finite-dimensional  approximations  of 
order  N  of  the  origineil  nonlinear  system  ^lre  considered  on  state  spaces  C  96,  and 
finite-dimensional  feedbacks  S^{x^^w)  are  computed.  These  feedback  functions  have 
the  same  structure  as  the  infinite-dimensional  functions  and  it  is  not  hcurd  to  write  down 
the  conditions  under  which  S{x,w)  as  N  -*  00.  It  is  not  enough,  however, 

that  all  of  the  individual  functions  in  the  feedback  formula,  e.g.  K  or  it*^,  converge.  In 
many  cases  we  must  require  that  these  functions  converge  in  a  locally  uniform  manner. 

Definition:  Let  (/„)  be  a  sequence  of  functions  defined  on  an  open  set  A.  We 
say  that  (/„)  converges  locally  uniformly  (l.u.)  to  a  function  f  if  every  a  ^  A  has  a 
neighborhood  J'f(a)  such  that  fn{x)-*f{x)  uniformly  for  x  in  J'f(a). 

The  condition  of  most  significance  for  controller  design  is  that  we  must  require  that  for 
every  ^€96,  G^(a)n^  •  ^  -»  G(a)-^  in  TL,  l.u.  in  a,  where  n^:96'^-»96  is  a  projection 
operator. 

To  proceed  further,  we  must  specify  how  the  state  feedback  gains  will  be 
computed  and  show  that  the  desired  convergence  does  indeed  occur.  In  this  work,  the 
feedback  gains  were  computed  by  solving  a  linear  quadratic  regulator  (LQR)  problem. 
Relying  heavily  on  the  previous  work  of  Gibson,  Banks,  and  Kunisch,  the  following 
result  was  proved. 

Let  96^,  jV  =  1,2,...  be  a  sequence  of  finite-dimensional  linear  subspaces  of  96 
and  n^:  96-^96^  be  the  canonical  orthogonal  projections.  Assume  that  for  every  a  €  -4., 
T^{t;a)  is  a  sequence  of  strongly  continuous  semigroups  on  96"^  with  infinitesim«J 
generators  A^{a)  €  i(96^).  Given  operators  B^{a)  €  i.(91, 96^)  and  D^{a)  G  X(96'^), 
we  consider  the  family  of  regulator  problems: 

Given  a  €  A.  minimize  J^{x^{0),u;a)  over  u  €  L^(0, 00;  *11),  where 

x^{t)  =  A^{a)x^{t)  +  B^(a)u(t),  t  >  0,  a  6  A, 

x«(0)  =  xS'  =  (1) 


and 


(2) 


Assumptions: 

HI  For  each  Xq  G  SB^and  for  each  a  €  -4.,  there  exists  an  admissible  control 
G  1/^(0, 00;*^)  for  (^Ri'^,a)  and  any  admissible  control  for  (1),  (2)  drives  the  state 
of  (1)  to  zero  asymptotically. 

H2  (i)  Every  a  G  -4  has  a  neighborhood  .^(a)  such  that  for  each  z  G  96,  we  have 
T^{t-,a)Il^ z—*T{t-,a)z  with  the  convergence  uniform  in  {t,cx)  on  I  xJ{{A),  where 
I  is  an  arbitrary  bounded  subset  of  [0, 00). 

(ii)  Every  a  G  -4.  has  a  neighborhood  >r(o)  such  that  for  each  z  G  96,  we  have 

z—*T{t;a)*z  with  the  convergence  uniform  in  {t,a)  on  7x>f(J.), 
where  I  is  an  arbitrary  bounded  subset  of  [0,oo). 

(iii)  For  each  u  G  'll,  {a)v—*B{a)v  locally  uniformly  in  a. 

(iv)  For  each  2  G  96,  B^ z—*B{a)* z  locally  uniformly  in  a. 

(v)  For  each  2  G  96,  z—*D{a)z  locally  uniformly  in  a. 

Theorem.  Suppose  (HI),  (H2)  hold,  Q{a)  >  0,  D{a)  >  0  and  D^{a)  >  0  and  let  R^{a) 
denote  the  unique  nonnegative  selfadjoint  Riccati  operators  on  96*^  for  the  problems 
(%^-,a).  Further  assume  that  for  each  a  €  A,  a  unique  nonnegative  selfadjoint  Riccati 
operator  on  96  for  the  problem  (%-,a)  exists  and  is  a  strongly  continuous  function  of  a. 
Let  S{t]a)  and  S^{t;a)  be  the  semigroups  generated  by  A(a)  —  B(a)Q~^(a)B*(a)i2(a) 
and  A^{a)  —  B^{a)Q~^{a)B^*{a)R^{a)  on  96  and  96^,  respectively  and  suppose 
II  5(t;Q)2  II  — »0  as  t—KX,  for  all  2  G  96  and  a  G -4.  If  for  every  a€A  there  is  a 
neighborhood  K{a)  and  positive  constants  M^,  M2  und  u  independent  of  N  and  t  such 
that 

\\S^{t-,a)\\^^<M^e-^^  fort>0,  N  =  l,2,...,  a  €  }{{a) 

and 

II  R^{a)  ^2  for  N  =  1,2,...,  a  G Jf (a) 

then 

R^{a)Il^z-*R{a)z  for  every  2  G  96,  locally  uniformly  in  a. 


This  result  replaces  the  conditions  on  the  LQR  feedback  gains  with  conditions  on 
the  open-loop  system.  To  verify  that  these  conditions  can  indeed  be  satisfied,  we 


restrict  attention  to  systems  whose  families  of  linearizations  consist  of  parabolic 
differential  equations  of  the  form 


i  =  a(s,  a)^-f  +  K-s*  +  c(s,  a)x  +  f{s,  a)u 


t>0,  0  <  s  <  1,  a  6  -4., 


with  initial  condition  x(0,  s)  =  boundary  conditions  r(t,0)  =  x(t,  1)  =  0,  and 

solutions  in  the  space  9G  =  1^(0, 1).  Define  the  operator  through  the  sesquilinear 

form 

<T{y,  z;  a)  =  I  [a(s,  -  (6(s,  +  {c(s,  a)  -  k)y)z  ]ds  =  {-  A^(a)y,  z) 

0 

where  A:  is  a  constant  chosen  to  malce  the  form  coercive.  Defining 

u,(a)  =  (/ - 

w*{a)  =  {I  -  Al{a))~^z, 

it  can  be  shown  that  H2  (i-ii)  is  satisfied  for  the  systems  being  considered  if  the 
following  condition  is  met: 


(Cl):  For  every  2€X^(0,1),  there  are  functions  w^:A-^'3o^  and 

N  =  l,2,...,  such  that  ||u;(a)-u;%)||  1-^  and  ||  u;*(a)  “  II  1-^  locally 

uniformly  in  a. 


H2  (iii-v)  are  easily  satisfied  by  appropriate  selection  of  the  approximating  problems. 
The  remaining  hypotheses  of  the  theorem,  save  the  strong  continuity  of  the  Riccati 
operator,  can  be  shown  to  follow  from  Cl  and  the  condition: 


(C2):  The  pair  (>l(a),  jB(a))  is  exponentially  stabilizable,  locally  uniformly  in  a.  This 
means  that  for  every  a  €  -4,  there  is  a  neighborhood  .N^(a)  and  a  bounded  linear  operator 
K s'.  independent  of  a,  such  that  for  every  a  €  ^(S),  the  semigroup  Ts{t\a), 

generated  by  A{a) -h  B{a)K g,  satisfies  IITg(t;a)ll  <  MjC  for  some  positive 
constants  Mj  and  Wj. 


To  show  that  the  above  theory  can  be  applied  to  a  relevant  example,  we  consider 
a  modified  form  of  Burgers’  equation  given  by 

Xi  =  eXgg  —  xxg  +  cx  +  uf 


and  subject  to  the  conditions 


x(5, 0)  =  Xo(,5),  x(0,  t)  =  x(l,  t)  =  Q,  0  <  5  <  1,  <  >  0. 

The  constant  c  is  used  to  make  the  equilibrium  point  at  the  origin  unstable,  u  is  the 
control  input  and  /  is  the  control  influence  function.  The  linearization  of  this  system  is 
given  by 

Xt  =  €Xss  -  x^{a)Xs  +  (c  -  x®(q!))x  +  a/, 

v.'here  x°{a)  is  the  equilibrium  state.  In  the  state  space  formalism,  we  identify  the  state 
operator  as 

A(q!)x  =  exss  —  x^{a)xg  +  (c  —  x2{a))x. 

The  linearizations  are  parabolic  partial  differential  equations,  as  required  by  our  results. 
Moreover,  { A(a);  a  €  ^}  constitute  an  analytic  family  of  operators  and  therefore 
generates  a  family  of  bounded  operators  T{t;a)  that  are  anailytic  in  t  and  a. 

The  finite-element  technique  is  used  to  obtain  finite-order  approximations  to  the 
original  system,  and  x(s)  is  approximated  as  a  piecewise-linecir  continuous  function. 
Relying  heavily  on  the  analyticity  of  T{t;a)  and  using  elementary  facts  concerning 
finite-element  approximations  we  have  shown  that  conditions  Cl  and  C2  are  satisfied  by 
our  example.  In  addition,  using  some  recent  results  of  Duncan  on  continuity  of  Riccati 
operators  depending  on  a  parameter  we  are  able  to  show  that  i2(a)  is  indeed  strongly 
continuous.  Thus,  we  can  be  assured  of  convergence  of  the  approximate  controllers  to 
the  ideal  controller.  A  number  of  computations  were  performed  to  illustrate  the 
convergence  numerically. 

In  summary,  we  have  placed  the  intuitive  idea  of  gciin- scheduling  on  a  rigorous 
mathematical  foundation  for  a  class  of  nonlinear  distributed-parameter  systems.  We 
have  considered  the  problem  of  approximating  the  ideal  controller  by  finite-dimensional 
controllers  cind  have  given  conditions  for  the  appropriate  convergence  to  take  place.  By 
specializing  to  the  case  of  ptirabolic  linearizations,  these  conditions  can  be  replaced  by  a 
much  simpler  set  of  conditions  which  we  have  shown  to  be  satisfied  by  a  modified  form 
of  Burgers’  equation. 

We  have  also  explored  several  interesting  problems  for  finite-dimensional 
nonlinear  systems.  Several  schemes  for  incorporating  global  information  into  a  gain- 
scheduled  controller  were  investigated.  As  discussed  above,  our  current  controllers  can 
not  be  guaranteed  to  produce  desired  performance  outside  a  neighborhood  of  the 
equilibrium  manifold.  To  modify  undesirable  performance  far  from  equilibrium,  it  is 
necessary  to  incorporate  the  desired  control  action  into  the  local  nonlinear  controller.  It 
Wtis  found  that  our  gain-scheduled  controller  provides  a  useful  structure  for 
incorporating  this  information.  Unstructured,  spline-based  methods  did  not  perform  as 


well.  We  have  also  explored  the  efficacy  of  a  nonliaear  PID  controller  applied  to  a 
highly  nonlinear  model  of  a  chemical  process.  The  performance  was  comparable  to  that 
achieved  by  much  more  sophisticated  nonlinear  controllers  that  required  more 
information  for  design  and  more  complex  on-line  computation. 
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